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The performance of the ILD detector at the ILC for the analysis of /i and r channels 
at the SUSY benchmark-point SPSla' has been studied with full detector simulation. 
It is concluded that if 500 fb"^ is delivered to the experiment, A{M^o) = 920 MeV/c^, 

A(A^Al) = 100 MeV/c^ A(M^o) = 1.38 GeV/c^ and A(cr(e+e- ^AlAl)) = 1-35 fb 
can be achieved from the /x channels alone. The preliminary results from the f channels, 
indicates that A(Pt) = 13 % is also achievable. 

1 Introduction 

The SUSY benchmark point SPSla' [2] offers a rich phenomenology at the ILC. It is point 
with quite low mass-spectrum in the slepton sector, and heavy squarks. Bosinos up to 
Xa (in e+e~ — ^XiXs) would be produced at Ecms = 500 GeV. It is a pure mSUGRA 
model, hence R- parity and CP is conserved. The unification scale parameters are: M1/2 = 
250 GeV, Mo = 70 GeV, = -300 GeV, tan (3 = 10, and sign{n) = +1. The point is not in 
contradiction with any experimental limits [5]. The f is the NLSP, and Mf^ = 107.9 GeV/c^ 
and M^o = 97.7 GeV/c^, so A(M) = 10.2 GeV/c^. At E^ms = 500 GeV, this yields 

min ^ GeV/c hence 77 events will pose a problem. As SPSla' is a point with 
an important co-annihilation contribution to the dark-matter relic density, the M^^ is a 
most important quantity to determine. An other consequence of the f being the NLSP, 
is that t:s are present in large fraction of the SUSY decays, so that SUSY itself will be a 
mayor background source for t channels. On the other hand, the Mfij^{Mfi,-^) is 189.9(125.3) 
GcV/c^, so that the minimum /i energy is 32.1(6.6) GeV. As, in addition, the branching 
ratios to jl in bosino decays are quite low, the fi final states offer cleaner conditions, and are 
well suited for doing the most precise measurements. 

The present note reports on the status of the analysis of of the ^ and r channels of 
the SPSla' scenario in the LDC detector. SPheno [4] was used to run the unification-scale 
model to the EW scale, and Whizard [5] was the used to generate events. The LDCPrime_02Sc 
detector model was fully simulated using MOKKA[B], and the events were reconstructed with 
MarlinReco [7]. The same chain was used to produce background events. 

2 Analysis of fj. channels 

Two channels containing only ^:s in the final state was chosen as a first study [5]: /ilAl 
HHXiXi ^nd X1X2 ~^ mArXi ~^ MMXiXi- As mentioned in the introduction, the SUSY 
background problem is not too severe in the /x channels, and it is advantageous to run 
the ILC at the polarisation giving the largest signal. Hence, these channels were studied 
assuming 80 % left e~ polarisation and 60 % right e+ polarisation. Under these conditions, 
the /UlAl process has a large cross-section, and is well suited to determine M^^ and M^a. 
X1X2 a small cross-section x BR, but can be used to determine M^a, without the need 
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Figure 1: The distributions of (left), and M^^ (right) 



to scan over the threshold. The main background processes are other SUSY giving two 
/i:s, mainly //r/xr, X2X2 with one X2 going to /i^, the other to i>u and ff with r ^iv^iVr 
Standard model background is mainly from WW and ZZ. Finally, each of the two processes 
is background to the other one. 

The following kinematic variables were used to disentangle signal and background, and 
to separate the two signal channels: The momentum of /i:s {P^), the acolinearity angle be- 
tween the ^:s (Oacoi), the acoplanarity angle between them, defined as the acolinearity in the 
projection perpendicular to the beam-axis {Oacop), the total missing transverse momentum 
(PTmiss) in the event, the invariant mass of the two fi:s (M^^), the total missing energy 
(Emiss), the polar angle of the missing momentum {Omissing p), and the velocity [3 of the 
H system. The distributions of P/j, and M^^ are shown in Fig. [TJ and that of /3 in Fig. [2] 

The fiLp-L channel was selected by de- 
manding that E^i,, e [200,430] GeV, 
M^^ ^ [80, 100] GeV and < 30 GeV/c^, and 
Omissing p £ [O.Itt, O.Ott] . Assuming an inte- 
grated luminosity of 500 fb~^, this leaves 
13000 events of SM background and 11000 
events of SUSY background, while 16300 
signal events were selected, corresponding 
to an efficiency of 60 %. The //l and Xi 
masses were then extracted by fitting the 
edges of the distribution, see Fig. [31 The 
errors on the fitted masses are A(Af^L) = 
100 MeV/c^ and A(M^o) = 920 MeV/c^, 
respectively. The beam-energy spread dom- 
inates these numbers. The production 
cross-section was determined using the extended likelihood formed by L{j)T^,9acoi))^ as 
these two variables were not used in selecting the signal. The uncertainty on the observed 
value is A((T(e"'"e~ -^/iLML)) = 1-35 fb. 

The X1X2 channel was selected by demanding that Omissing p G [0.27r, O.Stt], fi > 0.6, 
E„,,, e [355, 395] GeV, pxmiss > 40 GeV/c, M^^ e [40, 85] GeV/c^ and E„,, > 40 GeV/c^. 
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Figure 2: The distribution of the velocity f3 
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Figure 3: The fit to the lower (left) and upper (right) edges of the distribution 



At the assumed integrated luminosity, 
500 events of SM background and 2400 
events of SUSY background is expected, 
while 720 signal events were selected, cor- 
responding to an efficiency of 34 %. As- 
suming that the mass is known from the 
previous channel, the mass of X2 can be ex- 
tracted by a fit to the edge in the invariant 
mass spectrum, Fig. |4l and an uncertainty 
of A(M^o) = 1.38 GeV/c2 was found. 

3 T channels 

As mentioned in the introduction, SUSY it- 
self poses a background problem in the r 
analysis, and it is therefore needed to run the ILC at the polarisation that minimises the 
background. For 100 % left e~ polarisation and 100 % right e+ polarisation, the cross- 
sections for X2X2 ^"^cl XiXi s-re several hundred fb, and the branching ratios to f is above 
50 %. With the opposite polarisation, however, these cross-sections will almost vanish. 
Hence, these channels were studied assuming 80 % right e~ polarisation and 60 % left e+ 
polarisation. 

As the 77 background poses another challenge for the r channels, quite strong criteria 
must be applied: A correlated cut in p (the transverse momentum of the jets wrt. the thrust 
axis, in the projection perpendicular to the beam) was also done: p > 3 sin ^acop + 1.7. To 
further reduce the 77 background to acceptable levels, it was demanded that there be no 
significant activity in the BeamCal, and that the 4> angle of the missing momentum was not 
in the direction of the incoming beam-pipe. 

The f mass can be extracted from the end-point of the spectrum of Et- , which is equal to 
max' ^^'^ ^^^^ Xi mass, known eg. from the /iL analysis above. In principle, the maximum 
of the spectrum spectrum is at Pf so that the f can be used to find M^o as well, but 
due to the large 77 background, the maximum is quite hard to observe. 




40 50 60 70 80 90 



Invariant Mass [GeV] 

Figure 4: The fit to the M^^ distribution 
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Figure 5: E^et distribution after all cuts in the 
decay-mode independent analysis 



To extract the signal in order to de- 
termine the end-point the following cuts 
were appHed: Emiss £ [430,490] GeV, 
Mjet < 2 GeV/c^, Ojet above 20 degrees, 
Oacop < 160 degrees, 9acoi e [80, 170] de- 
grees, \cos9missmgp\ < 0.9, and charge of 
each jet — ±1. In addition, the anti-77 cut 
described above was applied. After these 
cuts, the SM background was 222 events, 
the SUSY background was 2747, while 8262 
signal events remained (10.2 % efficiency). 
Fig. [5] shows that the end-point is almost 
background free, and also that the turn- 
over point (expected to be at Pf = 
2.2 GeV/c) is too distorted by the cuts to be useful. 

The f mass-eigenstates are expected to be different from the chiral ones, and the off- 
diagonal term of mass-matrix is — Mt-(A^ — /itan/3). The diagonal terms in the mass matrix 
are known from ^-i^d M^j^ , so a measurement of O^ix gives Af ~ /i tan /3. If Xi is purely 
bino - it is in SPSla' - the r polarisation (Pt) depends only on 6mix- Ft can be extracted 
from spectrum for exclusive decay-mode(s). In this analysis, the r — > ir^ mode has 
been studied. The spectrum of 7r:s in the decay-chain f ^ r ^ tt^ Vt is shown in Fig. [SI 
with and without ISR and beam-spread. The highest sensitivity to the polarisation is in the 
region with P„ < P;^ 

The f ^ T ^ TT^ signal is selected 
with a set of cuts that intend to distort the 
spectrum as little as possible. The following 
pre-selections were first applied: The events 
should pass the anti-77 cut, Eyis should be 
< 120 GeV, the number of reconstructed 
particles < 20, and at least one of the two 
jets should contain a single particle. This 
single particle should be identified as a tt, 
and have E < A3 GeV. Finally, the to- 
tal charge should be 0. Events passing this 
preslection should then also fulfil the follow- 
ing criteria: The mass of the rest of the 
event (ie. after removing the signal pion) 
should be below 2.5 GcV/c^ | cos6\ of both 
jets should be < 0.9, and 6acop should be 

above 85 degrees. Finally, the sum over the two jets of the pT of one jet wrt. the direction 
of the other should be below 30 GeV/c. 

With these cuts, 134 SM jets remain, and 373 SUSY jets, while 2311 signal-jets are 
retained (13 %). The initial and final tt spectra are shown in Fig. [71 The procedure to 
extract the polarisation in the presence of background is to first fit the simulated back- 
ground alone to a heuristic function. The signal selections cuts are then applied to the 
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Figure 6: The P.^ distribution, at generator 
level. Shaded: fixed beam-energy, open: ISR 
and beam-spread included 



^When real data is available, the simulation of the background can be verified by reversing cuts to select 
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Figure 7: Distribution of before (left) and after (right) cuts. The right-hand plot also 
shows the spectrum after background-subtraction and efficiency correction (dots), and the 
final fit. 

signal-f background sample, and the function is subtracted from the observed distribution. 
An efficiency correction function, determined from signal-only simulation, is applied. The 
resulting distribution is then fitted with the theoretical spectrum, corrected for ISR and 
beam-spread, and the polarisation is obtained, see Fig. [71 right. Assuming an integrated 
luminosity of 500 fb^^, the value found is Pr = (93 ± 13) %, where the error also includes 
the uncertainty of the background parametrisation. 

4 Conclusions 

A study of some channels in SPSla' SUSY scenario fully simulated in the LDC detector 
at the ILC was presented. By analysing the channel e+c^ — >/iL/iL, it was concluded that 
A(M^o) = 920MeV/c2, A(M^J = lOOMeV/c^ and A(cr(e+e- ^mlMl)) = 1-35 fb, could 
be attained with an integrated luminosity of 500 fb^^ with 80 % left e^ polarisation and 60 
% right e+ polarisation. In the channel X1X2 ~^ f^P'RXi ~^ MMXiXii ^{^x"^ = 1-38 GeV/c^, 
was found, under the same conditions. It should be noted that this value is comparable to 
what a dedicated scan of the X2X2 threshold would give. 

In addition, a progress report on f production was given. The preliminary result on the 
measurement of the r polarisation gives A(P,-) = 13 %. Note, however, that this requires 
that the beam-polarisations are opposite to what the was used in the ^ channel. 

References 

[1] Presentation: 

http: //ilcagenda. linear collider . org/materlalDlsplay .py?contribId=173&sessionId=17&materialId=slldes&conf Id=2628 

[2] J. A. Aguilar-Saavedra & al., Eur.Phys. J.C46 (2006) 43, "arXiv:hep-ph/ 0511344l 

[3] C. F. Berger, J. S. Gainer, J. L. Hewett, T. G. Rizzo, arXiv:0812.0980 

[4] W. Porod, Comput. Phys. Commun. 153 (2003) 275, |arXiv:hep-ph/030l"l0T] 

[5] W. Kilian, T. Ohl, J. Router, arXiv0708.4233. 

[6] P. M. Freitas, MOKKA http : //mokka. In2p3 . f r/ 

[7] http: //ilcsof t .desy . de/portal/sof tware_packages/marlin/lndex_eng.html 

[8] N. d'Asccnzo, PhD thesis, DESY Thesis-2009-004 (2009) 

a signal-free, but SUSY-dominated region in the parameter-space. 



LCWS/ILC2008 



